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1. Introduction 

An accurate determination of the charm mass plays an important role on the precise physical 
evaluation of several observables, from K and B decays to CKM matrix elements and in lattice 
QCD. One of the usual techniques to extract the charm mass is to use the sum rules approach based 
on the relation between the moments of the production rate R and the inverse power of the square 
mass of the c quaik, and the Pade method (see ^ This approach should confront the fact 
that one have to employ the moments of the integral of R over the whole energy range, which are 
global properties, even though they are only known up to a certain scale A (since we only know 
experimentally 7? in a finite window). We propose to wield the local properties of R through a new 
"non-analytic reconstruction" method [^0. As we will show, this approach allows us to obtain 
local properties of the heavy quark correlators at each points of the spectrum with a systematic 
error and then to find a value of the charm mass directly on a regression on the experimental 
points. 



2. Details of the method 



2.1 Non-analytic reconstruction 

Let us consider the vector polarization function 

{qnqv-q^8nv) n{q^) = i J d'xe''''' (0|Tj^(^)/(0)|0) , (2.1) 

with the current j^{x) = \jf{x)Y^\i/{^)' which has a cut in the complex plane starting at = 4m^, 
where m is the (pole) mass of the heavy quark considered. In QCD perturbation theory, it can be 
expanded as 

where only IT^''' and IT'-^^ are know analytically, (for z = q^/4m^) 
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in which we used the auxiliary functions, 
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and 



(2.7) 



As it has been shown ^ even if the functions and n^^) are unknown analytically, one 
can reconstruct them from their expansions around <7^ — )• (Taylor expansion), — 4m^ (threshold 
expansion) and ^ oa (OPE), as 



NT 



n=Q 



n«(z) = £a«(«)«" + £(-)^' a«LiW(p,«)-i8^5LiW +^«(iV,*,«). (2.8) 



Let emphasize a little this expression. First one defines the so-called conformal change of 
variable 



4(0 
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(2.9) 



This change of variables maps the cut z plane into a unit disc in the CO plane, as we can see on 



Figure 2A_. The physical cut z G [l,oo[ is transformed into the circle |g)| = 1 . The points z = into 
ft) = 0, z = 1 into ft) = 1 , the limit z — > +°° ± ie into ft) — > — 1 it /£, and z — > —°° into ft) — ^ — 1 . 
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Figure 1: Conformal mapping between z and (O. 



For both functions n^^' and n'^', Feynman diagrams calculations at — )• give the expan- 
sions up to an order Nl (for k = 2,3) 
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where the relation between the two coefficients C^*^' and H^*^' {n) is 

c«(«) = 2i-2«r(2«) y — ^- 



(2.11) 
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The main part of the approximation in (2.8) Ues on the combination of the poly logarithms 
functions, 



LiW(5,w) 
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and the analytic evaluation of the coefficients \ and jS^ \. In order to reconstruct n'^' and n'^), 
we collect here their corresponding coefficients (see [Rl H|| for more details) 
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At least, one gives the error functions S'^'^\ 
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which encode the systematic error from the reconstructions . 

2.2 Experimental data 

There exists several experimental results for the e^e^ in hadrons that one can use for the fitting 
of the c quark mass. Each of the experiments give the ratio R{s) of the radiation-corrected measured 
hadronic cross section to the calculated lowest-order cross section for muon pair production, 



R(s) 
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(2.18) 



that has the experimental values shown in Fig. ^ . 

This Fig ^ shows that the complete spectrum is sensitive to resonances, as expected. It is 
obvious that a perturbative approach cannot take into account the resonances description, then one 
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Experiment 


Reference 


MARK I 
PLUTO 
CrystalBall (Run 1) 
CrystalBall (Run 2) 
MDl 
CLEO 
CLEO 
BES 
BES 
CLEO 
CLEO 


© 
[0] 
§] 
ii 
t 

[12] 
[13] 
[14] 
[15] 



Table 1: All different experimental sets considered for the fits. 



MARKI 1986 ^ 
PLUTO — 
ChiystalBall 1990 ^ 
MD-1 1996 ^ 




5 10 15 20 25 

(in GeV) 



Figure 2: Collection of the different experimental sets for the V-V spectrum. 

has to make an arbitrary choice on where we assume that the continuum limit is reached or in 
other words, where the perturbative description is pertinent. Let's choose the value of 5 GeV. Of 
course the influence of the arbitrariness has to be discussed and taken account in the evaluation of 
the error but it is something depending on the perturbative and heavy-quark limit more than the 
reconstruction itself. 

The idea now is to peiform a fit among all this data points to extract the perturbative mass 
of the c-quark. 
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2.3 Fitting approach 

The first step in the fitting procedure is to choose the following expression for the running 
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where A is the energy scale and the /3 -function has coefficients 
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and is the number of quarks with mass smaller than y/s/2. 
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The theoretical expression (2.18) is related to Tl{q^) ( [2.2[ ), up to a^, 

a, (5) 
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where all n^*^) functions have the argument z = and Nc is the number of colors. 

The goal of the analysis is to extract from the comparison between the value of /?exp. and 
/?th.- The usual method used is to built the moments associated to R from to and identifying 
the coefficients of the Taylor expansion that are proportional up to a factor to m^^. Instead of this 
approach, we propose to perform the analysis directly on the function itself, because thanks to the 



reconstruction method formula (p.8[), its expression is available and its systematic error too ( |2. 16[ ). 
For this we will use a ;^^-method with the assumption 



(2.22) 



where the sj are the experimental energy points, the a^xp. is the experimental error and the theoret- 
ical error ath. due the approximation of the reconstruction is given by 
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3. Results 

3.1 Numerical results at order aj 

At aj order, one obtains after a regression procedure with a ;^j^;j^/d.o.f. = 1.03, 

nicipole) = 1.85 ±0.08 GeV , 



that is translated into the MS mass as [ ]16| ] 

mc(MS) = 1.12 ±0.08 GeV 



(3.1) 



(3.2) 



;^min/d-O.f=1.03 

mdpole) = 1.85 ± 0.08 £FiV 



Assuming now that the mass obeys to a Gaussian density of probability, one can easily 
reconstruct points by points the eiTor generated on /?th. by this hypothesis, taking into account that 
the relation between and /?th. is highly non linear and non trivial for expressing the eixor. We 
choose then to use a Monte-Carlo approach to obtaining the mean value of 7?th. and its error as 
shown in Fig ||. 
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Figure 3: The reconstructed radiation-corrected measured hadronic cross section to the calculated lowest- 
order cross section for moon pair production. 



4. Conclusions 

We show that it is possible to extract the charm mass value after a regression to the exper- 
imental data of the radiation-con^ected measured hadronic cross section to the calculated lowest- 
order cross section for muon pair production using the non-analytic reconstruction of the heavy- 
quark correlators. We present here a preliminary result up to a^. The next step would include the 
order a| and a complete analysis of all different systematic contributions [17|. 
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